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PURPOSE 


The  firm  of  Lee  Arnold  Associates  presented  a proposal  to  the  Air 
Force  which  woiild  initiate  an  experimental  investigation  of  the  merits 
of  a novel  helicopter  rotor  system.  It  was  recommended  at  that  time 
that  a theoretical  investigation  he  conducted  in  order  to  justify  the 
initiation  of  such  an  experimental  program.  The  results  of  the 
theoretical  study  are  contained  herein. 


CONCLUSIONS 


It  is  concluded  from  the  foregoing  that  the  use  of  the  Kdrmin-Yuan 
rotor  incorporating  an  oval  shaped  airfoil  section,  boundary  layer 
suction,  and  flap  control  on  a helicopter,  would  make  possible  substan- 
tial improvements  in  performance  characteristics,  including  irate  of 
vertical  climb,  forward  flight  speed  and  useful  load  lifting  capacity. 

It  is  further  concluded  that  an  experimental  investigation  is 
warranted  in  order  to  substantiate  the  findings  presented  in  this 
report  and  to  obtain  Information  that  could  be  applied  to  an  actual 
design  incorporating  the  Karman-Yuan  rotor  features. 
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SUMMm 


A performance  calculation  has  been  carried  out  for  the  purpose  of 
compeurlng  the  Plaseckl  XE>l6  rotor  vlth  the  Kim^Ln-Tuan  rotor.  The  tvo 
rotors  vere  assumed  Identical  In  size  and  lifting  capacity  (rotor 
diameter  « 82  ft,  gross  vel^t  » 38600  lbs).  The  Plaseckl  rotor  blade 
Is  equipped  vlth  an  VkCk  OOI5  alxfoll  and  the  Kirmin-lban  rotor  blade 
vlth  an  elliptical  airfoil  of  l8^  thickness  chord  ratio  emd  cyclic  flap 
control.  The  controllable  lift  distribution  along  the  blade  span  of  the 
Karman-Yuan  rotor  Is  achieved  by  varying  the  spanvlse  flap  location  and 
applied  boundary  layer  suction.  Due  to  the  high  lift  coefficient 
enqployed  In  the  K&rm^-Yuan  rotor  blade  It  vas  found  that  a lower  rota- 
tional speed  can  be  used  to  sustain  the  same  gross  velfi^t  as  does  the 
XH-16. 


The  results  of  the  calculations  are  summarized  In  the  folloving 
tables . 


Table  1 Plaseckl  XH-I6  Helicopter  (tip  speed  63O  ft/sec) 


Rotor  Horsepower 

Induced 

Profile 

Parasite 

Total 

Hovering 

1984 

740 

0 

2726 

160  MPH 

410 

843 

1415 

2668 

Table  2 KAmAn-Yuan  Helicopter  (tip  speed  430  ft/sec) 


Rotor  Horsepower 

Induced 

Profile 

Parasite 

Suction 

Total 

Hovering 

2020 

137 

0 

250-350 

2407-2507 

160  MPH 

410 

179 

1234 

250-350 

2073-2173 
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The  rotor  suction  horsepower  includes  an  equivalent  power  loss  due 
to  weight  of  suction  equijanent.  The  range  of  values  of  that  term  res\ilts 
from  assuming  maximum  and  minimum  estimates  corresponding  to  laminar  and 
possible  turbulent  flow  on  the  rotor  blade. 

It  is  seen  from  the  above  that  the  saving  in  rotor  horsepower  for 
the  hovering  case  of  the  Karman-Yvian  rotor  (over  that  of  the  Piasecki 
XH-16  rotor)  is  about  IO56  and  the  power  saving  at  16O  MPH  is  about 
These  favorable  characteristics  result  from  the  lower  rotational  speed 
eii5>loyed  and  from  the  lower  profile  drag  coefficient  on  the  blades 
achieved  by  suction  of  the  boundary  layer.  This  lOjt  saving  in  hovering 
flight  corresponds  to  a 305^  increase  in  payload  for  the  helicopter  under 
consideration.  Alternatively,  for  the  same  gross  weight  the  rate  of 
climb  of  the  Karm^n-Yuan  helicopter  would  exceed  the  Piasecki  XH-I6 
by  8036. 

In  addition  to  bringing  about  a large  power  saving  within  the  for- 
ward flight  range  of  the  Piasecki  XH-I6  helicqpter,  the  K^rm&n-Yuan  rotor 
would  eliminate  the  present  speed  limitation  of  the  conventional  heli- 
copter which  would  result  from 

(1)  stalling  of  the  retreating  blades,  and 

(2)  compressibility  effects  in  the  advancing  blades. 

The  effect  described  in  the  above  item  (2)  is  due  to  the  low  tip  speed 
and  low  peak  of  the  potential  velocity  on  the  upper  blade  surfaces. 

An  aeroelastic  investigation  indicated  that  certain  unusual  aero- 
elastic  characteristics  were  found  to  be  inherent  in  the  Thwaites  airfoil, 
which  it  is  believed  could  be  circxanvented  by  suitable  structural  design. 

I 

It  should  be  noted  that  extremely  conservative  estimates  were  made 
in  the  computation  of  iKJwer  requirements  of  the  K^rm^n-Yvian  rotor.  It 
is  felt  that  somewhat  less  conservative  estimates  would  reveal  sharply 
improved  operating  characteristics. 


WADC  TR  53-497 


viii 


IWTRODUCTION 


V 


1- 


This  report  presents  the  results  of  an  Investigation  including  a 
performance  study  and  an  aeroelastic  investigation  of  a Piaseckl  XH-I6 
helicopter  equipped  with  a KaLrman-Tuan  rotor.  The  performance  study 
included  the  effect  of  replacing  the  conventional  rotor  system  of  the 
XH-16  helicopter  by  a Kirm&n>Yuan  rotor. 

The  XH-I6  is  a twin  rotored  helicopter  having  82  foot  diameter 
rotors  arranged  in  overlapping  tandem  emd  its  physical  dimensions  and 
aerodynamic  design  characteristics  are  as  follows: 


Ifumber  of  blades 
Blade  chord  (constant) 
Blade  radius 
Rotor  solidity 
Blade  airfoil  section 
Gross  weight 
Tip  speed 
Power  available 


3 per  rotor 
2.33  feet 
41  feet 
.0546 
MCA  0015 
36654  pounds 
636  feet/second 
3020  horsepower 


The  Kerman- Y\ian  rotor  is  provided  with  a doubly  synmetrlc  blade 
section  with  a cyclic  flap  control  (Patent  Serial  No.  3^6,326,  1953)* 

The  periodic  variation  of  lift  on  the  blade  in  forward  flight  is  ful- 
filled by  cyclic  variation  of  flap  position.  The  production  of  high  lift 
is  then  Independent  of  blade  angle  of  attack  and  the  separation  of  the 
flow  is  prevented  by  continuous  suction  of  the  bOTindary  layer.  The 
physical  dimensions  of  the  helicopter  with  Karman- Yuan  modifications  are 
kept  the  same  as  the  XH-16  but  the  difference  in  aerodynamic  design 
characteristics  resulting  fr<mi  use  of  higher  lift  coefficients  are  as 
follows: 


Blade  airfoil  section  18^  elliptical  section 

Tip  speed  430  feet/ second 


In  the  vertical  flight  condition  the  average  lift  coefficient  for 
the  XH-I6  rotor  blade  is  .442  suid  for  the  Kaxman-Yuan  rotor  blade  it  is 
assumed  to  be  1.0.  This  is  feasible  because  of  its  non-stall  character- 
istics. Due  to  the  higher  lift  coefficient  eoployed  in  the  Kirmdn-Yuan 
rotor  blade  a 32.51^  decrease  in  rotational  speed  can  be  used  to  sustain 
the  same  gross  weight  as  the  XE-I6.  In  the  comqtarison  of  vertical 
flight  perfozmemce  it  is  sufficiently  accurate  to  assume  that  the  rotor 
blades  in  both  cases  are  Ideally  twisted. 
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It  is  shown  in  Section  5 that  the  induced  powers  consumed  in  any 
two  rotors  (identical  in  size)  are  approximately  inversely  proportional 
to  their  respective  tip  loss  factors.  Since  the  difference  in  tip  loss 
of  the  above  considered  rotors  is  less  than  2$  their  difference  in 
Induced  power  is  very  slight. 

A greater  saving  in  profile  power  consvunption  by  the  Karman-Yuan 
rotor  over  that  of  the  XH-I6  results  from  the  lower  rotational  speed 
employed  by  the  former  rotor  and  from  its  lower  profile  drag  coefficient 
on  the  blades  achieved  by  suction  of  the  boundary  layer.  The  Justifica- 
tion of  the  profile  drag  coefficient  used  in  the  calculation  of  profile 
power  for  the  KArmAn-Yuan  rotor  was  explained  in  Section  5*  Because  of 
this  greater  saving  in  profile  power  of  the  Kaiman-Yuan  rotor  over  that 
of  the  XH-I6  a great  Increase  in  vertical  rate  of  climb  is  achieved  by 
the  former  helicopter  system. 

In  the  forward  flight  performance  calculation  especially  at  the 
cruising  speed  range  the  Induced  power  consumption  is  only  about  2^^  of 
the  total  power  consun^jtion.  Because  of  a greater  saving  in  profile 
power  as  well  as  parasite  power  due  to  the  rotor  blades  by  the  Karman- 
Yuan  rotor  over  that  of  the  XH-I6,  the  total  power  consumption  in  for- 
ward flight  is  much  less  than  that  of  the  XH-I6.  The  method  used  in 
forward  flight  performance  calculation  does  not  take  into  account  the 
loss  in  the  region  of  reverse  flow.  Since  the  reverse  flow  has  no 
detrimental  affect  on  the  KArmAn-Yuan  rotor,  it  is  evident  that  the  per- 
formance coa^Mirlson  in  the  forward  flight  condition  is  conservative. 

The  suction  power  was  calcxilated  so  as  to  be  sufficient  to  compen- 
sate for  all  losses  occixrring  in  the  flow  system.  In  order  to  take  into 
consideration  the  possibility  of  the  flow  becoming  tviAulent  the  actual 
suction  quantity  used  in  the  present  calculation  was  arbitrarily  assumed 
to  be  ^ times  higher  than  the  calculated  values.  Extremely  conservative 
assun^tlons  were  made  throughout  in  estimating  suction  requirements. 
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a 
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Circulatory  force  (Equation  (8.33)) 

Suction  pressure  due  to  lift 
Duct  pressure  drop  due  to  loss  of  momentum 
Pressvire  drop  in  duct  due  to  friction 
Duct  pressure  due  to  centrifugal  force 
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Equation  (8.26) 

Profile  torque  or  suction  quantity 

Induced  rotor  torque 
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Suction  quantity  through  porous  skin  per  unit  radial  length 
of  blade 

Badlus  of  rotor  disc 
Effective  radius  of  rotor  disc 
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(Figure  8.2) 
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2 

(Vhov)  A^c  (eq^'v^lent  downwash) 

Total  weight  of  helicopter 

r/R  or  coordinate  of  a point  in  direction  of  major  axis  of 
airfoil 

Effective  working  blade  range  (Equation  (5 *3)) 

2 

(e  Uo)/2aV 
Angle  of  attack 

Flap  angle  (subscript  t for  tip  angle) 

Rotor  tilt 


xlv 


WADC  TR  53-^97 


r 

A 

f 

r 

*Ponc 

'fonc 

e 

To 

Jl 

JL 

uO 


Measuz^  of  boundary  layer  thickness 

Displacenent  thickness 

(VoX  )/0 

(®/v)  an/ds 

Density  of  alrstream 

Rotor  solidity 

Inflow  angle  of  blade  element  (subscript  t for  tip  angle) 

Hon-clrculatory  potential  on  contour  of  airfoil 

Hon-clrculatory  stream  fvinctlon  on  contour  of  airfoil 

Momentvim  thickness  or  polar  angle  of  point  on  contoiir  In 
^ - plane  (Page  41) 

Wall  shearing  stress 

Angular  velocity  of  rotor  (Sections  5,  6,  7)  or  of  blade 
section  (Section  8) 

Vector  angular  velocity  (Equation  (8.17)) 

Frequency  of  oscillation  of  blade  section  (rad. /sec.) 
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SECnOH  1 


SURVEY  OF  UTBRATURE 


1.1  OlLCrr  TESTS 


A report.  Reference  (1),  on  a wind  tunnel  test  of  a 35^  thick 
elliptical  airfoil  has  been  studied.  The  high  lifting  characteristics 
of  the  Thmltes  type  airfoil  vere  substantiated,  for  a naxlnun  lift 
coefficient  of  7*3  ws  obtained. 

However,  this  Investigation  was  very  lijslted;  no  Beasurements  of 
drag  vere  made,  the  suction  was  applied  over  the  total  surface  of  the 
airfoil,  and  was,  consequently,  greatly  In  excess  of  requirements. 
Consequently,  no  estimate  of  optimum  suction  power  requlrmnents  was  made 
available.  Effects  of  flap  chord  variation  and  chord  location  vere  not 
Investigated. 


1.2  BRITISH  IHVESTIGATIOHS 


The  British  reports  (References  (2)  emd  (3)  ) contain  the  under- 
lying principles  of  the  Thvaltes  type  airfoil.  This  Includes  a method 
for  detexmlnlng  profiles  which  could  obtain  xmlform  chordvlse  lift 
distributions. 

However,  there  Is  practically  no  experimental  Information  available 
except  for  some  two-dimensional  tests  on  circular  cylinders,  the  purpose 
of  idilch  was  to  substantiate  the  basic  principles.  It  was  seen  that  the 
high  lift  coefficients  predicted  by  theory  could  be  realized  and  that 
wake  drag  was  practically  eliminated. 

Ho  suction  meae\irements  vere  made. 
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1.3  UHITED  AIRCRAFT  REPORTS 


Various  United  Aircraft  Reports  vere  studied  (References  (4),  (5), 
(6)  and  (7))*  Only  References  (4)  and  (5)  contained  pertinent  Irfonna- 
tlon. 


Sunmarles  of  these  reports  follow: 


1.  U.A.C.  Report  No.  R - 278,  Reference  (4).  This  report  presents 
an  investigation  of  the  optimum  arrangement  for  maximum  lift  and 
Jet  thrust  of  various  slot  and  flap  combinations.  The  results 
of  tests  of  a two-dimensional  MCA  23012  wing  model  equipped 
with  pressure  slots  can  be  sisnmarized  as  follows: 


l.a  A maximum  lift  coefficient  of  4.7  was  reached  with  a ^ 

chord  plain  slot  at  75^  chord,  a ^ chord  lip  type  slot  at 
100^  chord,  flap  angle  of  6o°,  and  wing  pressure  of 
l6x(  j?  )r/2) . 

l.b  An  Increase  of  quarter-chord  moment  coefficient,  from 
C]b1.£  -0.26^  for  the  basic  wing  with  flap  aud  no  slots 
(Cl  may  - 2.14)  to  Cjji  = -1.10  for  wing  with  flaps  and 
slots  (Cl  tuny  » 4.7)  IS  indicated. 


l.c  The  suction  quantity  coefficient  Cq  s .0094  was  obtained 
Cl  may  = 3*72,  and  Cq  « .03I  was  obtained  for 
Cl  msuc  ■ ^•7* 


l.d  It  is  pointed  out  that  the  Increased  mcmentum  due  to  the 
Jet  may  reduce  the  drag. 


2.  U.A.C.  Report  No.  R - 95374-3#  Reference  (5).  This  report 
presents  an  e:q>erlfflental  investigation  of  the  effect  of  trail- 
ing edge  suction  and  pressure  slots  to  the  lift  control  of  a 
wing.  The  purpose  of  this  investigation  is  to  apply  this  lift 
control  scheme  as  a gust  alleviation  device.  The  results  of 
the  investigation  can  be  sinmnarlzed  as  follows: 

2. a For  the  purpose  of  alleviating  the  effect  of  most  gusts 
for  cm  average  3C0  nph  flight  a lift  coefficient  of  0.4 
would  be  sufficient. 

2.b  The  most  efficient  configuration  tested  was  the  trailing 
edge  pressure  slot  which  produces  a lift  coefficient 


r 


A 
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increment  of  about  0.4  for  a power  coefficient  of  .005. 
Such  a power  coefficient  would  represent  about  7 to  of 
the  total  airplane  power.  However,  in  order  to  produce  a 
lift  coefficient  increment  of  about  1.4,  a suction  power 
coefficient  of  0.6  is  needed. 

2.C  The  pitching  mcmient  increments  were  neeurly  Independent  of 
all  variables  except  lift  Increments.  The  value  for 
d Cm/d  Cl  is  about  -O.3  which  is  very  close  to  the  value 
obtained  in  the  same  case  of  pressure  slots  and  flaps. 


It  is  believed  that  the  mechanisms  described  in  the  above  reports 
would  not  be  suitable  for  application  to  helicopter  rotors  because  of 
the  large  moment  coefficients  as  ccnpared  .to  other  hi£^  lift  devices  and 
high  suction  quantity  coefficients  accoopEUiylng  the  increase  of  lift 
coefficients . 


1.4  CONCLUSIOHS 


The  reports  studied  do  not  fiumish  the  information  that  is  required 
to  evaluate  the  characteristics  of  the  Thwaltes  type  airfoil,  as  applied 
to  helicopters,  there  are  no  other  data  available  to  acccnpllsh  this 
evalmtlon. 
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SECTION  2 

VELOCITY  DISTRIBUTION  ABOUT  AN  ELLIPTICAL  AIRFOIL 


The  velocity  distribution  for  the  potential  flow  around  an  ellipti- 
cal airfoil  can  be  calculated  frcm  the  known  complex  potential  for  a 
circle  with  aid  of  the  following  transfoxnatlon  between  a circle  and  aui 
ellipse  (see  Reference  (8)). 

? z * ) 2-1 


where  c^  - a^  - b^  and  the  expression  z ^ cosh  ^ represents  elliptical 
coordinates^*,  x = a cos/^  and  y = b sin  ^ (on  the  ellipse). 

If  the  potential  flow  consists  of  a parallel  and  a circulatory  flow, 
the  velocity  at  any  point  on  the  elliptical  airfoil,  as  determined  by 
the  above  mentioned  transformation  can  be  written 

U 0+  + 


The  value  of  s/2a  correspondlog  to  a given  value  of 
from  the  Integral 


Is  determined 


which  Is  easily  evaluated  by  means  of  elliptic  Integrals. 


The  coordinates  of  the  stagnation  points  for  any  given  C];,  amd  o( 
can  be  detezmlned  by  the  following  expressions: 


2^ 

<x 


cos 


C Cl 

2Tr(l+^)  * 4-Tr^(l+  ^/af' 


1 

h 


- Cl.  COt,oC 

2Tr(‘+  ^/a) 


QCj 


Cl 

4Tr'^(i+ 


ZA 


The  nomenclatiires  In  Equations  (2.1),  (2.2),  (2.3)  and  (2.4)  are  ex- 
plained In  Figures  2.1  and  2.2. 

The  velocity  distributions  for  the  potential  flow  aroxmd  an  elllptl 
cal  airfoil  of  17*36^  thickness  chord  ratio  and  Cl  =•  1>  2,  4 and  c<  = 0 
are  calculated  and  shown  In  Figures  2.3>  2.4  and  2.5,  respectively.  The 
corresponding  first  derivative  of  the  velocity  distributions  are  given 
In  Figures  2.6,  2.7  and  2.8,  respectively.  It  will  be  noted  that  the 
velocities  and  gradients  are  given  only  for  the  forward  half  of  the 
ellipse  (top  8uid  bottom)  since  they  are  respectively  symmetric  and  anti- 
symmetric with  respect  to  the  vertlcsd  center  line. 
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SECTION  3 


CALCUIATION  OF  SUCTION  QUANTITY 


In  the  previous  section  a potential  flov  is  calculated  about  an 
elliptical  airfoil  for  any  given  lift  coefficient  and  the  location  on 
the  airfoil  of  the  dividing  streamlines  is  determined.  If  a flap  is 
placed  along  the  calculated  rear  dividing  streamline  (at  the  stagnation 
point),  the  steady  flow  of  real  fluid  past  the  airfoil  will  approximate 
closely  to  the  calculated  potential  flow  with  circulation  (References 
(1)  and  (2)),  In  order  to  prevent  separation  which  would  destroy  the 
high  lift,  continuous  suction  must  be  applied  along  the  surface  of  the 
airfoil  where  a positive  pressure  gradient  exists. 

If  a one-parameter  family  for  the  velocity  profile  in  the  boundary 
layer  is  chosen,  the  following  differential  equation  for  non-dimensional 
momentum  thickness  is  obtained  from  the  Khrman  Momentum  Equation 


■ (^) 


[a-  2Ag(k) 3. 


where 


2 • 

0"  Uo 

'Za'? 

3(k)- 

® i 

^ ^ + 0. Co  ^ 

II 

(t  -2)  k 

A - "A 

K - 

1-^(1-  f 

A = 

s'  ^ 

V <^s 

V 

- O.OZ5  5 8 K' 


3.Z 
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The  differential  equation  for  momentim  thickness  without  suction 
can  be  obtained  from  Equation  (3.1)  by  putting  /^  = 0 and  K = A -1, 
which  gives 


2g(A)  f Ca) 


3.1, 


where 

g(A)  = = 0.403S  -H  O.II57A-  0.02358A^ 

3. 

f Ca)  = 0-523<<>  - /.4S3  A - o*''2'73  a’" +■  .04-7  I6>  A^ 


From  the  expressions  for  A and  Z , (Equation  (3.1))  and  (Equation 
(3.2)),  we  obtain 


A f (A) 


3.  S’ 


Equations  (3*3)  and  (3.5)  are  the  parametric  representation  of  a 
first  order  differential  equation  in  Z( s ) . They  may  be  solved  by  pro- 
cediires  such  as  the  isocline  method  provided  the  velocity  distribution 
U and  dU/ds  and  the  initial  values  of  Ao  and  ( dz/ds  )o  at  the  forward 
stagnation  point  Sq,  are  known.  It  is  shown  in  Reference  (8)  that 


3.  (p 
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The  numerical  solution  of  the  pair  of  Equations  (3»3)  smd  (3*5)  can 
be  considerably  sliiq)lified  by  eliminating  the  parameter  A • lo&y 

done  by  plotting  2fg  versus  A 8^  fitting  straight  line  segments  to 
various  portions  of  the  curve.  It  was  found  that  the  single  straight  line 


2-9  (J  = ,47  - 6.53  A 3*^ 


3,7 


is  sufficiently  accurate  in  the  range  of  A encountered.  Substituting 
Equations  (3.3)  and  (3*5)  into  (3*7)  w©  have  the  first  order,  linear 
differential  equation 


yz 

6.53 

This  has  the  solution 


3.9 


The  second  term  in  the  bracket  is  zero  when  Sq  is  a steignatlon  point, 
but  must  be  Included  in  other  non- suction  regions. 

The  solution.  Equation  (3*9)  is  obtained  by  simple  quadratures  in 
the  region  of  positive  velocity  gradient,  du/ds  > 0.  At  the  mlnlmm 
pressura  point,  at  tdilch  the  potential  velocity  gradient  becomes  zero, 
dU/ds  = 0,  the  velocity  profile  inside  the  boundary  layer  approaches  the 
flat  plate  velocity  profile  (K  = -1).  If  the  parameter  K in  Equations 
(3*1)  and  (3»2)  is  kept  constant  emd  equal  to  -1,  the  velocity  distri- 
bution in  the  boxmdary  layer  will  remain  constant  along  the  body  in  the 
flow.  Equation  (3‘1)  can  then  be  reduced  to  the  form 

3.10 

it) 


which  has  the  slnq)le  integral 
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where  Si  is  the  distance  along  the  body  to  the  Initial  point  at  which 
the  suction  is  begrin  and  Zi  is  the  corresponding  value  of  Z,  determined 
from  the  solution  of  Equation  (3*9) • 

The  above  calculations  are,  of  course,  valid  only  iq)  to  the  transi- 
tion region,  as  Equation  (3*1)  does  not  apply  to  turbulent  flow.  In 
regard  to  the  transition  Reynolds  Humber  at  which  flow  may  become 
tiurbulent,  it  can  only  be  said  at  this  time,  that  due  to  the  favorable 
velocity  profile  secured  by  use  of  suction,  the  critical  Reynolds 
Dumber  becomes  much  higher  than  for  the  case  without  suction.  The  above 
interpretation  is  based  on  the  theory  of  the  stability  of  the  laminar 
boundary  layer,  as  in  Reference  (10). 

Equation  (3*9)  gives  the  non-dimensional  momentum  thickness  as  a 
function  of  velocity  in  the  potential  flow  and  the  initial  values  at 
^ich  the  application  of  suction  begins.  Since  all  the  items  are  given 
in  non-dimensional  form,  the  value  of  the  non-dimensional  momentum 
thickness  calculated  from  Equation  (3*9)  is  unrestricted  as  to  free 
stream  velocity,  length  of  chord,  and  the  viscosity  of  the  air.. 


For  a given  thickness  to  chord  ratio  of  an  airfoil,  the  non-dimen- 
sional velocity  distribution  in  the  potential  flow  can  be  c(mq?uted  as  a 
function  of  the  non-dimensional  arc  length.  The  non-dimensional  momentum 
thickness  in  the  region  vhere  suction  is  not  needed  is  calculated  from 
Equation  (3*9) • The  suction  is  then  applied  at  the  minimum  pressure 
point,  and  Equation  (3*11)  is  used  to  calculate  the  mcmientum  thickness. 
The  initial  values  Zj,  (U^/Uq)  and  (8/2a.)i  in  Equation  (3«11)  occur  at 
the  end  point  of  the  solution.  Equation  (3'9)> 


The  suction  velocity  Vq  necessary  to  maintain  a constant  shape  of 
the  velocity  profile  in  the  boundary  layer  can  be  derived  from  Equation 
(3.2),  which  gives 


0 


l<4t 

3.IZ 
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where  R,  is  the  Reynolds  Number  with  reference  to  the  chord  of  the  air- 
foil. The  value  for  9 in  Equation  (3-12)  is  obtained  from 

Equation  (3*9) • 

The  resulting  curves  of  momentxim  thickness  versus  the  arc  length  of 
the  airfoil  are  shown  in  Figures  3*1^  3*2,  3 *3;  and  the  suction  velo- 
cities are  shown  in  Figures  3*4,  3*5»  3*6  for  2,  4,  respectively. 

The  suction  quantity  needed  in  order  to  maintain  the  flat-plate 
velocity  profile  in  the  boundary  layer  after  the  minimum  pressure  pointy, 
can  be  detemined  from  Equation  (3*12).  The  coefficient  of  suction 
quantity  Cq  can  then  be  fovind  by  the  following  expression 


Equation  (3.13)  can  be  expressed  in  the  form 


A maximum  value  of  kj.=  3.33  will  be  used  in  Section  7.2  to  obtain  the 
suction  quantity  Qq.  This  value  os  obtained  from  Figures  3.4,  3.5, 
and  3*6. 
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SECTION  4 


CALCULATION  OF  PROFILE  DRAG 


The  surface  profile  drag  on  the  given  airfoil  can  be  calculated 
from  the  solution  of  laminar  boundary  layer  theory  given  in  the  previous 
section.  From  the  following  momentum  equation  for  the  botmdary  layer 
with  pressure  gradient  and  suction 


V 


die 


i! 

e 


(2+A)eu 


dO 


= uv. 


r 


4.1 


we  have 


U dyf? 


A-,1 

For  Vq  ^ 0 and  K = -1,  Equation  (4.2)  is  reduced  to  the  following 
expression 


— _ o>42<?2.  fo  \ 


4.-5 


For  Vq  = 0,  the  following  equation  for  friction  drag  is  obtained. 


^(Vo)  Il^AL- 
v/T 


4.4 


The  coefficient  of  friction  drag  is 


1 

^ Up(>er 

‘iu.f'feii  e. 

LouJC*' 

WADC  TR 

53-497 

11 

4S 


Subfitituting  the  results  of  Equations  (4.3)  and  (4.4)  into  Equation  (4.5) 
yields  the  following  numerical  results. 


Cl 

1 

2 

4 

Cdo 

iM 

4^ 

6.00 

rwi 

■TrI 

The  above  results  are  plotted  In  Figuire  4.1.  A linear  relationship 
Is  indicated  between  the  profile  drag  and  lift  coefficients,  and  a value 
of  = 2.55  is  thereby  found  for  Cl  - 0. 

It  must  be  borne  In  mind  that  the  above  estimates  are  based  on  leun- 
Inar  flow  theory.  It  Is  likely  that  part  of  the  flow  in  the  neighborhood 
of  the  rear  stagnation  point  might  become  tvirbulent  and  that  the  actual 
friction  coefficient  may  become  two  to  three  times  greater  than  the 
above  calculated  values. 

In  order  to  justify  the  above  conclusions,  the  seune  method  was  used 
to  calculate  the  profile  drag  coefficient  for  the  HACA  0015  airfoil.. 

For  laminar  flow  with  suction  It  gives  j^Cd^  = 3*02  for  Cr  =0.5.  The 
elliptical  airfoil  considered  in  the  present  report  gives  =3*0 

for  the  same  lift  coefficient. 

The  profile  drag  coefficient  of  the  NACA  0015  airfoil  was  then  cal- 
culated for  the  case  of  zero  suction  and  consequent  turbulent  boimdary 
layer  in  the  region  of  positive  pressure  gradient  with  the  result 
Cdo  ~ .0107  for  Rj  = 4.7x10°.  For  a Reynolds  Number  very  close  to  that 
value,  the  experimental  value  of  the  profile  drag  coefficient  for 
HACA  0015  is  .00928,  thereby  showing  close  agreement. 

In  view  of  the  above  comparison  it  is  reasonable  to  assume  that  the 
calculated  values  of  the  profile  drag  coefficients  may  be  multiplied  by 
a factor  of  3*5  in  the  estimation  of  performance  of  the  KArmAn-Yuan 
rotor  when  the  elllptlccJ.  airfoil  is  used  in  order  to  teUce  into  account 
the  possibility  of  turbulent  flow. 
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SECTION  5 

STUDY  OF  THE  EFFECT  OF  THE  KARMEN- YUWf  ROTOR  MODIFICATION 
ON  TEE  PIASECKI  XH-I6  HELICOPTER  VERTICAL  FLIGHT  PERFORMANCE 


In  the  present  section  a comparison  of  vertical  fli^t  performance 
between  the  Plaseckl  XH-I6  helicopter  and  the  proposed  Karman-Yuan 
modification  was  made.  Since  there  is  no  test  data  available  on  the 
Karm&n-Yuan  rotor,  it  is  not  Justified  at  the  present  to  make  a detailed 
compsorison  such  as  effects  at  various  altitudes,  and  Mach  Number  correct- 
ions. It  is  sufficient  to  assume  that  the  blades  are  Ideally  twisted 
for  both  above  mentioned  rotors,  and  that  they  are  both  operated  at  sea 
level  conditions . 


5.1  RELATION  BETWEEN  THRUST  AND  INDUCED  POWER 


The  differential  thrust  expression  for  a rotor  having  b blades  can 
be  obtained  from  the  blade-element  theory  as  follows. 

o(T  - dr  5./ 

For  the  case  of  ideal  twist 

Cu--  ^ 'ft) 

where  ^ is  the  flap  angle  defined  in  Section  2 and  ^ is  the  Inflow 
angle  of  a blade  element.  The  subscript  ( )-t  represents  the  condition 
at  the  tip  of  the  blade,  c and  a are  respectively  the  chord  and  dCL/doc  . 

Substituting  Equation  (5.2)  into  Equation  (5.1)  and  replacing  r/R 
by  X,  the  total  expression  for  thrust  becomes 
Xe 
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The  limits  of  the  integral  are  taken  vithln  the  effective  working  blade 
range  to  Xg.  Equating  Equation  (5.3)  to 

T-  C-rirR^fC-aR)’' 


the  expression  for  thrust  coefficient  becomes 
Ct  - -?«■(. 

If  Cl  represents  the  average  lift  coefficient  on  a blade  the  expression 
for  thrust  coefficient  may  be  written  as  follows. 

Ct  - f C,  S.S 

The  inflow  angle  at  the  tip  of  a blade  can  be  obtained  from  the 
combined  momentum  and  blade-element  theories  which  gives 


S.(f> 


The  flap  angle  at  the  tip  of  a blade  is  obtained  from  the  Equations 
(5.4)  and  (5*6)  as  follows. 


4-Ct  / Ct 

(ra(xj-x;)  2Cx^-x>) 


5.7 


The  induced  part  of  the  rotor  torque  which  is  caused  by  the  compo- 
nents of  the  lift  vectors  in  the  plane  of  rotation  is  written  as 


o 

l^on  substitution  of  and  from  Equations  (5.4)  and  (5.6) 

into  Equation  (5-5)  the  coefficient  of  Induced  torque  becomes 


a ( Xg-  X- )"/»' 


5.  *7 
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The  induced  power  per  rotor  is  given  by  the  following  expression. 


P 


I 

550 


C«;  j>TrR'(_aRf 


5JO 


5«1.1  Induced  Power  Calculation  for  Piasecki  XH-I6  Rotor 
5.1.1  (a)  Rotor  characteristics 


Gross  weight 
Tip  speed 
Blade  chord 
Rotor  diameter 
Rotor  solidity 


36654  (twin  rotors) 
636  ft/ sec 
2.33  ft 

82  ft 

.0546 


5.1.1  (b)  The  thrust,  torque  and  Induced  power  coefficients 
are  computed; 


W/z 


.003<p2. 


Xe=  « - 


\f£o^ 


‘111 


,Z 


- _ ^ I 
" <r  CX|-X?) 


,44Z 

.000((<>l(p 


The  Induced  power  of  a single  XH-I6  rotor  (7^6 
blade  twist)  was  calculated  in  Reference  (l2)  to 
be  1004  HP.  The  variation  of  the  blade  angle  and 
the  inflow  angle  along  the  blade  are  shown  in 
Figure  5.1, 
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5.1.2  Induced  Power  Calculation  for  KAimAn- Yuan  Rotor 

In  the  present  conqparlson  the  average  lift  coefficient  of 
1.6  is  assumed  for  the  elliptical-flap  airfoil  system.  With 
this  average  lift  coefficient,  the  thrust  coefficient  for 
the  rotor  is 

^ = .00^01  x|  - .OOOOJZCp 

Substituting  = 4.5(  1 - ) into  the  above  equation  a 

cubic  equation  for  Xg  yields  the  following  solution. 

Ct  = .007^1 

In  order  to  keep  the  same  rotor  thrust  and  rotor  dieuneter, 
the  tip  speed  of  the  Kiimin-Yuan  rotor  becomes 


_aR 


T 

CrpTTR^ 


43  0 ft/^ec. 


as  compared  to  JIR  = 636  ft/sec  for  the  XH-16, 
The  induced  torque  coefficient  is 


C 


r x’’ 

Ct  ^e 


.O0O5B75 


and  the  induced  power  per  rotor  is 

^ = 1012  hP 

The  variation  of  the  flap  emgle  and  the  inflow  angle  along 
the  blade  are  shown  in  Figvire  5*2. 

It  will  be  noted  that  the  ratio  of  the  Induced  power  of  two 
different  rotors  can  be  expressed  by  the  following  simple 
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expression 


P.2 


-Xi) 


Approximately,  the  Induced  power  ratio  of  two  different 
rotors  Is  Inversely  proportional  to  the  ratio  of  their  tip 
loss  factors.  Since  the  difference  In  the  tip  loss  ratios 
Is  less  than  2^,  the  difference  In  Induced  power  Is  also 
less  than  2^. 


5.2  PROFILE  POWER 


The  differential  profile  drag  expression  for  a rotor  having  b 
blades  Is 

dO^  ^ bf  S-'/ 

where  Is  the  profile  drag  coefficient  of  the  blade  section. 

Thus  the  differential  power  required  to  overcome  the  profile  drag 
of  b blade  elements  Is 

dPo  = SJZ 


5.2.1  Profile  Power  for  XH-16  Rotor 

In  Reference  (12)  a value  of  Cdo  ~ .00928  was  used  to  cal- 
culate the  profile  power  of  the  XH-l6  rotor  giving 

= 3*7^  I 
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5.2.2  Profile  Power  for  the  Kann^-Yuan  Rotor 


The  profile  drag  coefficient  of  em  elliptical  airfoil-flap 
system  vas  calculated  in  Section  4 and  it  caji  be  expressed 
by  the  following  equation 


^Jo  ~ 
where 

c,  -- 


2.55  t 0.^53  Cu 

C^t?)c  _ C43o)(233) 

V ‘ lOx,o-‘  ^ 


CS.I3) 


Since  the  profile  torque  is  given  by 


or 


Qo- 


(T 

T~ 


I 

irj)R’(-aR)*j'Cj„X*c(X 


C5.I4) 


the  profile  torque  coefficient  becomes 


.00/5  Y 


cs-is9 
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Hence  the  rotor  profile  drag  coefficient  of  the  KarmAn-Yuan 
rotor  Ifi  .0015. 

In  order  to  be  conservative  the  profile  drag  coefficient 
used  in  the  calculation  of  the  profile  pover  vas  assumed  to 
be  .0055  thereby  accounting  for  possible  turbulence. 

We  obtain 


Po-  ^ ^ a ^.S  iP/r.ior 


5.3  SUCTIOM  POWER  FOR  THE  kArM^N-YUA-N  ROTOR 


The  power  required  for  boundary  layer  suction,  sufficient  to  main- 
tain a given  flow  condition  and  lift  distribution  along  the  rotor  blades 
is  discussed  in  detail  in  Section  7*  Tbe  total  suction  power  required 
per  rotor  is  found  to  be  175  HP. 


5.4  SUMMARY  OF  ROTOR  HORSEPOWER  FOR  HOVERIHG 


Pi 

Po 

Psuction 

Ptotal 

Reference 

Section 

(5.1) 

(5.2) 

(7.2) 

Flaseckl 

XH-16 

992 

371 

0 

1363 

KirmAn-Yusm 

1012 

68.5 

175 

1255.5 

it  of  Power  Saving  by  Kiro^-Yuan  Rotor  _ 1363-1255 .5  - 9»2?6 

1363 
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The  above  power  saving  is  equivalent  to  a increase  in  payload 


5.3  COMPARISON  OF  VERTICAL  CLIMB  PERFORMANCE 


The  vertical  rate  of  climb  is  computed  by  comparing  the  total  flow 
through  the  disc,  U,  with  the  downwash  velocity  v.  The  difference 
between  these  two  values  is  the  rate  of  climb. 


It  is  shown  in  Reference  (12)  that  the  rate  of  vertical  climb,  Vg, 
is  related  to  the  induced  downwash  in  hovering,  vj^Qy,  and  to  the  total 
flow  through  the  disc  in  vertical  climb,  U^,  by  the  approximate  formula 


-L 


5.16 


where 


U 


c 


550 


^AV  ~ 

w 


and 


55  O 


Pi 

W 


Therefore  the  vertical  rate  of  climb  in  feet  per  minute  becomes 


Vt  = 60  (it  - 


where 


1. 

"^hov 


5.17 


The  resulting  characteristics  of  the  two  rotors  are  listed  below. 
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XH-16  Helicopter  (Twin  Rotors) 


Pi 

Po 


^hov 


^av 


U. 


1984  HP  itrom  Section  6.1.1) 

742  HP  (from  Section  6.2.1) 

550  Pi/W  = 29.8  ft/sec 

3020  HP  (from  Reference  (12) , available  HP) 

550  (Pav  - Po)  = 34.2  ft/sec 
W 

(^ov^)/Uc  = 26  ft/sec 
60  (Xfc-  Vgq)  = 492  ft/min 


Kara^-Yuan  Helicopter  Modification  (Twin  Rotors) 

Pi  2024  HP  (from  Section  6,1*2) 

Po  137  HP  (from  Section  6.2.2) 

Pav  3020  - 350  = 2670  HP 

VhoY  550  Pi/W  = 30.4  ft/sec 

"^eq  ('^hov^)/Uc  - 24.3  ft/sec 

Vc  60  (U^-  Vgq)  = 821  ft/min 

^ increase  in  Vc  by  Karman-Yuan  Helicopter  Modification  1=  19*5$ 
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5,6  POWER-OFF  DESCEMT 


In  case  of  power  failure,  it  is  necessary  for  a helicopter  to  glide 
in  autorotation  with  a low  descent  speed.  In  general,  the  blade  angle 
of  attack  as  well  as  the  blade  lift  coefficient  required  for  autorota- 
tion must  be  con^ratively  small.  This  can  be  accomplished  in  the 
Karman-Yuan  rotor  by  decreasing  the  angle  ^ to  a small  magnitude 
i.e.,  by  moving  the  flaps  close  to  the  trailing  edge  of  the  blades. 

Since  the  lift  coefficient  is  coiiparatively  small,  the  suction  required 
can  be  relatively  small.  It  appears  that  this  small  suction  can  be 
accomplished  as  a result  of  the  combined  centrifugal  and  tip  venturi 
effects  discussed  in  Section  J.  Since  the  rate  of  descent  of  a heli- 
copter depends  to  a great  extent  on  the  rotor  profile  drag,  it  is  clear 
that  the  lower  profile  drag  of  the  Karman-Yuan  rotor  resulting  from 
suction  would  lower  the  rate  of  descent. 
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SECTION  6 


STUDY  OF  THE  EFFECT  OF  THE  KARMAN-YUAN  ROTOR  MODIFICATION 
ON  THE  EEASECKI  XH-I6  HELICOPTER  FORWARD  FLIGHT  PERFORMANCE 


The  method  used  in  the  calculation  of  forward  flight  performance  is 
based  on  the  outline  given  in  References  (12)  and  (l^).  Since  this 
method  does  not  take  into  account  the  loss  of  efficiency  in  the  reverse 
flow  region,  and  the  K^rmAn-Yuan  rotor  does  not  suffer  this  loss,  it  is 
evident  that  results  calculated  in  the  present  section  are  conservative. 


6.1  INDUCED  POWER 


The  induced  power  in  the  region  of  forward  speeds  of  particular 
interest  (from  that  of  best  climb  to  Vnax)  represents  only  a relatively 
small  fraction  (1/3  to  1/5)  of  the  total  power  required,  hence  the 
method  based  on  momentxam  theory  is  sufficiently  accurate  for  practical 
prediction  of  the  induced  power. 

The  Induced  power  which  is  required  to  compensate  for  the  energy 
carried  away  by  the  slipstream  when  downwash  is  produced  can  be  expressed 
in  horsepower  as 


S50 


I 


where  W is  the  weight  supported  by  the  rotor  and  v is  the  induced 
velocity.  By  analogy  with  fixed  wing  aircraft  the  average  downwash  may 
be  found  as  a function  of  the  area  of  the  cross  section  of  a stream  tube 
having  as  its  diameter  the  rotor  span  2R.  Thus,  for  a single  rotor  heli 
copter  the  average  downwash  velocity  is 


tT  = 


W 

ZTTR’-pV' 


<p.2 


where  V is  the  resultant  flow  through  the  disc  and  may  be  expressed  as 


V'  X + nT  + 
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where  ^ is  the  tilt  of  the  rotor  disc.  However,  since  for  the  region 
of  low  ftwd  moderate  flying  speeds  is  small  in  congiarison  with  v, 

while  for  high  flying  speeds  the  influence  of  all  axial  coii5)onents 
(v  and  ) are  mimportsait  compared  to  Vj^,  for  all  practical,  purposes 

the  resultant  V can  be  expressed  approximately  as 


V'  - 6.3 

When  the  rotors  are  arranged  in  close  tandem,  however,  it  is 
evident  that  the  average  induced  velocity  will  increase  due  to  their 
proximity.  The  thrust  of  a tamdern  rotored  helicopter  may  be  eaqpressed 

0>.4- 

where  Ae  is  the  effective  projected  area  of  both  rotors,  is  the 
effective  rotor  radius,  and  Vjj  is  the  horizontal  component  of  the  total 
air  flow  and  Gap  is  the  vertical  distance  between  rotor  disc  planes. 

For  forward  speeds  above  60  nph  Eqmtlon  (6.4)  can  be  simplified  by 
letting  , which  results  in 


The  resiilts  of  the  calculation  of  induced  power  for  the  XH-I6  heli- 
copter and  the  Karman-Yuan  modification  are  listed  below.  (Both  heli- 
copters have  approximately  the  same  Pi*) 


Re  = 41  ft. 


A^:  ‘?<p<5>2  J “ 3.5<e  -ft 


/iT  = - 24-. G \/^, 


33.24- 


<o.S<K 


From  Equation  (6,1)  we  obtain  the  induced  power 
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(mph) 

Vf, 

(fps) 

V 

(fps) 

(®) 

60 

88 

15.18 

1007 

80 

117 

11.88 

788 

100 

147 

9.77 

649 

120 

176 

8.42 

559 

140 

206 

7.07 

470 

160 

234 

6.17 

410 

6.2  PROFILE  POWER 


In  Reference  (13), it  bae  been  shown  that  the  profile  power  in  for- 
ward flight  can  be  closely  approximated  by  increasing  the  hovering  pro- 
file power  as  a function  of  the  square  of  the  advance  ratio  M.  , Thus 
the  profile  power  in  forward  flight  becomes  ' 


(i*/) 


Since  the  rotor-tilt 
practice. 


angle  ^ is  usually  small 


yL 

-X2.R 


is  used  in 


The  results  of  the  calculations  of  the  profile  power  for  the  XH-I6 
and  the  Karman-Yxian  modification  are  listed  below. 


\ (ii5>h) 

XH-16 

KarmAn-Yuan 

n 

1^/ 

Po  (HP) 

e 

1.0 

742 

1.0 

137 

40 

1.0085 

748 

1.0186 

140 

60 

1.0192 

756 

1.0418 

143 

80 

1.040 

772 

1.0746 

147 

100 

1.0532 

782 

1.1162 

153 

120 

1.0766 

799 

1.170 

160 

140 

1.1041 

820 

1.232 

169 

160 

1.1359 

843 

1.302 

179 
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6,3  PARASITE  POWER 


The  total  parasite  drag  of  the  whole  helicopter  is  composed  of  the 
parasite  drag  of  the  helicopter  body  and  that  of  the  rotor.  The  rotor 
contribution  to  the  parasite  drag  may  be  esipressed  as  follows 


I 

zw 


be  ( JI  r + 


ip  j^sin  \p 


" pev  r-o+or 


(p7 


The  profile  drag  coefficient  of  a conventional  rotor  is  higher 
than  would  result  from  the  drag  analysis  of  the  individual  blades  because 
of  the  occurrence  of  reverse  flow  areas.  We  will  here  also  assume  a 
slightly  higher  to  cover  the  effect  of  possible  flow  spoilage  in 
reverse  flow,  although  that  may  not  be  required. 

The  parasite  drag  of  the  helicopter  body  (except  rotors)  will  be 

where  the  equivalent  flat  plate  area  f of  the  XH-I6  helicopter  body  is 
as  given  in  Reference  (12). 

The  total  parasite  power  of  a tandem  rotored  helicopter  is  express- 
ed as  follows : 

"iT3o 


6.3.1  Calculation  of  Parasite  Power  for  XH-I6  Helicopter 

f = 39.4  ft^  0“».0546  R ^ 41  ft 

- .014  J^R  -636  ft/sec 

« 2.16x10-6v2(39.4v^+-  2570) 

V h n 

The  resTilts  are  stmsnarlzed  in  the  following  table. 
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6.3 « 2 Calculation  of  Parasite  Power  for  Itiirman-Ymn  Helicopter 
Modification 

f =39.4  ft^  r*.0546  R = 4l  ft 

Cdo  = .0085  I2.R-  430  ft/i 

Pp=  2.16x10"^V^^(39.4V^+  1055) 


(mph) 

^PB 

(HP) 

^PR 

(HP) 

^p 

(HP) 

40 

17.2 

8.0 

25 

6o 

58.0 

17.7 

76 

80 

137.5 

31.3 

169 

100 

269 

49 

318 

120 

464 

71 

535 

140 

735 

96 

831 

l6o 

1108 

126 

1234 
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6.4  SUCTIOH  POWER  FOR  KABM^-YU(UT  HELICOPTER  MODIFICATION 


The  eviction  power  needed  In  forward  flight  is  actually  lees  than 
that  required  in  hovering  as  a result  of  the  small  lift  coefficient 
needed  in  forward  flight.  However,  the  suction  power  used  below  is 
assvmed  to  be  the  same  as  in  hovering,  i.e.,  350EP  for  both  rotors. 

A summary  of  forward  flight  performance  of  the  two  systems  is  listed 
below. 


(mph) 

■ 

40 

60 

80 

100 

120 

140 

160 

\ 

(fj)8) 

0 

58.7 

88 

117 

147 

176 

206 

234 

Pi 

Induced 

1984 

1335 

1007 

788 

649 

??? 

470 

410 

PO 

Profile 

742 

748 

756 

772 

782 

199 

820 

843 

Plaseckl 

XH-16 

mm 

0 

36 

101 

214 

388 

636 

968 

1419 

pt 

Total 

m 

n 

1864 

1774 

1819 

1^ 

2258 

2668 

Pi 

Induced 

— 

2020 

_.1335 

1007 

788 

649 

??? 

470 

410 

Po 

Profile 

137 

140 

143 

147 

153. 

160 

169 

- 179 

Karman- 

Yuan 

wmm 

0 

25 

76 

169 

n 

m 

HI 

mi 

Ps 

Suction 

350 

350 

350 

350 

350 

350 

350 

..  350 

Pt 

Total 

1850 

1576 

1454 

1470 

1604 

1820 

2173 

Power  Ratio 
Karm^-Yuan/XH*l6 

.92 

^ 

^85 

.82 

.81 

00 

• 

.81 

.81 

Curves  of  total  power  versus  speed  are  plotted  in  Figure  6.1. 
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SECTION  7 


CALCUIATION  OF  SUCTION  POWER  REQUIREMENTS 


7.1  EXPRESSION  FOR  SUCTION  PRESSURE 


A suction  pump,  or  an  ejector  and  venturi  tube  at  the  blade  tip,  is 
required  in  the  contemplated  Karman-Yuan  rotor  to  provide  sufficient 
suction  pressure  to  draw  a desired  quantity  of  air  mass  against  the  max- 
imum localized  negative  pressure  increments  over  the  vlng  surface.  In 
addition,  the  suction  provided  by  the  p\nq>  must  also  be  sufficient  to 
overcme  the  pressure  drop  due  to  duct  friction,  loss  of  momentvan,  and 
the  pressure  drop  across  the  porous  skin.  In  this  section  the  term 
suction  pressure  will  be  used  to  denote  pressure  difference  between  the 
atmospheric  pressure  and  the  absolute  pressure,  p^  - p. 

In  a uniformly  rotating  system  such  as  a helicopter  rotor  an  addi- 
tional component  of  force  in  the  direction,  of  flow  within  the  duct  which 
is  due  to  the  component  of  the  centrifvigal  acceleration  SlP'  r in  that 
direction  must  be  added  to  the  above  mentioned  system.  The  coriolis  . 
acceleration  is  always  at  right  angles  to  the  direction  of  the  relative 
flow  and  hence  yields  no  contribution. 

The  equation  for  fluid  motion  in  rectangular  ducts  can  be  written 
as  follows: 


where 

U^  s duct  velocity  at  section  r, 
p^  - suction  pressure  at  section  r, 
p^  r suction  pressure  due  to  lift, 
h = the  depth  of  rectangular  section  duct. 
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b,  - the  width  of  rectangular  section  duct,  and 
= duct  wall  shearing  stress. 

Integrating  Equation  (7«1)  ve  obtain  an  expression  for  the  duct  suction 
pressure  at  any  section  r 

Pj  = p*  * p,  *pi- 

where 


and 

fL  = 

4 2(b  + h)  hydraulic  diameter) 


Since  the  volume  flow  in  the  duct  is  U^b^h  and  the  volume  flow  into 
the  duct  through  the  porous  surface  in  an  element  of  length  dr  is  . 

Vgb,  dr,  the  following  continuity  equation  can  be  established. 

^ (Uj  b,h)  r -jj 

For  a duct  with  constant  cross-section  the  duct  velocity  becomes 

Uj  = 7,4 

where  Qg  is  the  suction  qmntity  per  \inlt  radial  length. 


P3 

H-  X2 


SO 
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From  Section  3 we  obtain  the  following  relation 


s- 


and 


-five 

V 


But 


Q<,  ~ )c  = J _fz.rc9 


7.5 


where  kj.  is  a function  of  the  lift  coefficient  Cj^.  Substituting  Equation 
(7.4)  into  Equation  (7*3)  yields 


/-TLCT? 


Jk,fr 


dr 


If  is  assumed  to  be  the  maximum  value,  say  k,  then 


d 3 f\ 


> ^ - J 


^ere 


^ J^-CLR}cV  j Jt  - length  of  duct 
From  Equation  (7«7)^  the  exJ)reBsion  for  beccMes 

The  expression  for  P2  becomes  ^ 

/Qo\*  •*  * / ^ \ 


T#  <p 


7.7 


7.r 


7.^ 
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vhere 


Cf  = .00(4  + .OIZ5 

J^c 

and 


/ 


l?c 


SLPfC. 

V 


The  eai)re8sion  for  Pj^  Is  computed  using  the  maximum  value  of  U at  each 
section  in  order  to  be  conservative. 


7./0 


•»diere  Pg  is  the  pressxire  across  the  porous  skin.  The  expression  for 
P3 1. 


-Q.  r 


T.M 


7.2  SUCTION  DUCTS  AITO  SUCTION  QUANTITY 


Since  the  distribution  of  the  lift  pressure  and  suction  quantity 
varies  in  the  chordwlse  direction  it  is  more  economical  to  use  several 
suction  ducts  than  a single  one.  However,  in  the  present  estimate  only 
a single  duct  is  considered  which  is  assumed  to  take  all  the  suction 
quantities  at  the  most  adverse  condition  of  wing  pressure  loading. 

The  area  of  the  cross  section  of  the  duct  is  assumed  to  be  one- 
third  of  the  blade  cross-section  area  (6.256  ft^)  and  it  is  located  in 
the  rear  part  of  the  blade.  The  calculated  values  of  the  perimeter  of 
the  duct  is  2.235  ft.  Thus  the  hydraulic  diameter  is 
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The  suctloDnquantlty  Is  obtained  f rom  Cq  (page  lO); 

Q^=  3.33  V = 1.333 

therefore 

Oo  _ C 'o 
"A  ■ 


In  the  present  calculation  a factor  of  2.5  Is  multiplied  to  the 
calculated  value  of  Qq/a 


With  the  above  parameters  the  various  pressures  \dilch  contribute  to 
the  duct  pressures  can  be  calculated  as  follows: 


0002.1  +.0011^ 

A-  2 ^■2- 

l<g.T3  X + (.774  X 


3.SZ 


Jp-i  " ^ - 220  x'*' 

where  u/_a.r  versus  x is  shown  in  Figure  7.1.  The  distribution  of  C.  is 
obtained  by  considering  the  linear  variation  of  flap  angle  & versus  x so 
that  we  can  also  plot  Cl  versus  U/.fl.r  in  Figure  7.1.  ' 
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The  computations  are  s\iBmiarized  in  the  following  tables. 


(x) 

(x)2 

(x)3 

(x)3-52 

i-(uiir)2 

0.2 

.04 

.008 

.00343 

.0016 

- 9.72 

0.3 

.09 

.027 

.0144 

.0081 

- 5.68 

0.4 

.16 

.064 

.0395 

.0256 

- 3.5 

0.5 

.25 

.125 

.0895 

.0625 

- 2.33 

0.6 

.36 

.216 

.165 

.1296 

- 1.56 

0.8 

.64 

.512 

.456 

.410 

- .99 

0.96 

.922 

.885 

.866 

.850 

- .69 

(x) 

Px 

■^3 

^d 

0.2 

85.2 

2.4 

- 8.8 

78.8 

0.3 

112.5 

8.1 

- 19.8 

100.8 

0.4 

123 

19.2 

- 35.2 

107 

0.5 

128 

37.5 

- 55 

110.5 

0.6 

124 

64.8 

- 79.2 

109.6 

0.8 

125 

154 

- I4i 

138 

0.96 

140 

265 

- 203 

202 

Pg  is  negligible 


The  calculated  mEiximum  duct  suction  pressure  at  the  tip  is 
20^/ft^.  The  designed  maximum  suction  pressiure  in  the  duct  is  assumed 
to  be  300  ^/ft^  in  order  to  be  consenrative.  The  required  suction 
power  for  one  blade  is 


3k 
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p. 


(3oo)(4/) 

S50 


Q Hdx. 

5 

• “Jo 


= 4-2  IP 


; ^5  ^ Qo  \r^ 


(30  0)  (4/) 
S5^0 


(l3 


-9^ 


.-2. 


Total  suction  for  one  rotor  is  126  HP,  (S^-Ps  three  blades), 

%rt  “ t " |(l^)r+l>(.'»3«)  = 331  «/ssc. 

^9*  ,S6. 

<?T  = CAOCl^x.^-Stp)-!- [x  = 8S  cu.^f^ec 


If  the  weight  of  the  blower  is  estimated  at  4 poxmds  per  horsepower 
and  the  duct  equipment  is  assumed  to  be  100  po\mds  per  rotor  the  total 
weight  of  the  suction  system  is  about  600  pounds  per  rotor.  For  the 
XH-16  helicopter  the  power  loading  is  about  13  pounds  per  horsepower. 
There  will  be  required  an  additional  power  of  46  horsepower  per  rotor 
due  to  the  weight  of  suction  equipment.  Therefore  the  total  suction 
power  required  is  about  175  horsepower  per  rotor  for  the  maximum 
estimate. 

This  estimate  can  be  reduced  considerably  if  a tip  venturi  or  tip 
ejector  is  used  to  draw  the  boundary  layer  in  the  outer  part  of  the 
blades  and  a pump  system  is  used  to  draw  the  boundary  layer  in  the  inner 
portion  of  the  blades. 
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7.3  REVERSE  FLOW  CHARACTERISTICS 


One  of  the  limiting  factors  in  performance  of  orthodox  helicopters 
is  that  resulting  from  the  loss  of  lift  in  the  region  of  reverse  flow 
when  a helicopter  is  in  forward  flight.  It  is  believed  that  this  loss 
of  lift  and  limiting  effect  on  performance  can  be  avoided  to  a great 
extent  by  application  of  the  Karman-Yuan  system. 

If  an  additional  flap  is  placed  along  the  dividing  streamline  in 
the  vicinity  of  the  leading  edge  of  the  oval  shaped  airfoil,  approxi- 
mately the  same  lift  would  be  produced  as  there  would  be  without  the 
leading  edge  flap.  Now,  if  the  flow  over  the  airfoil  is  reversed, 
approximately  the  same  lift  would  again  be  obtained,  for,  aside  from 
separation,  it  can  be  shown  that  the  pattern  of  streamlines  and 
pressure  distribution  is  independent  of  the  sign  of  the  velocity  of  flow. 
The  only  requirement  is  that  in  addition  to  the  suction  at  the  trailing 
edge,  suction  must  be  applied  at  the  leading  edge  in  that  region  where 
reverse  flow  would  be  experienced.  For  this  reason,  the  Karman-Yxian 
rotor  blades  wo\ild  not  experience  loss  of  lift  in  the  reverse  flow 
region.  Furthermore,  the  higher  blade  lift  coefficient  obtained  at  the 
inner  portion  of  the  blade  would  result  in  a favorable  spanwlse  lift 
distribution. 
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In  the  analysis  of  the  divergence  of  a rotor  blade,  the  torsional 
resisting  moment  (r)  developed  by  the  blade  as  a restilt  of  a rotation 
, is  a rather  con^dicated  function  of  the  radial  distance  r along  the 
blade.  Nevertheless  ve  may  define  a section  torsional  stiffness  Co((r) 
such  that 

= <C^Cr) 

Equating  the  total  aerodynamic  moment  at  section  r with  the  resisting 
moment  , we  obtain 


or 


1 

Cc 

II 

23f  Sinf 

20C  ^ 

(o.-  yyi  J 

If  the  dencmlnator  of  Equation  (8.4)  Is  zero  or  negative  for  every 
section  of  the  rotor  blade,  then  the  divergence  velocity  for  the  entire 
blade  must  become  Infinite  or  Imaginary.  The  divergence  criterion  can 
be  expressed 


5in  2< 
Z«. 


Za.  (o(+e) 

<x  - b «=<. 


^ O 


or 


q-b  sin  < _ - 
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j’rcm  Equations  (8.1)  and  (8.2)  it  is  seen  that  the  right  hand  side 
of  the  inequality  (8.5)  is  the  location  of  the  center  of  pressvire  a^^,  so 
that  the  divergence  criterion  (8.5)  requires  the  elastic  axis  to  be 
located  upstream  of  the  center  of  pressure  in  order  to  preclude  the 
possibility  of  divergence  at  any  speed. 

If  it  is  assximed  that  is  a small  variation  from  the  initially 
zero  angle  of  attack,  the  Equation  (8.5)  can  be  vritten 


a 


g-b  oc‘ 

2 sm  I 


We  may  conclude  that  for  zero  initial  angle  of  attack,  (center  of 
pressure  at  the  center  of  the  airfoil),  location  of  the  elastic  axis  at 
or  upstream  of  the  center  of  the  airfoil  insures  an  infinite  divergence 
speed. 

The  above  criterion  for  divergence,  vhich  involves  the  location  of 
the  center  of  pressiure  of  the  total  lift  on  the  airfoil  is  inadequate  in 
discussing  stability  of  the  airfoil  with  regard  to  perturbation  lift 
forces.  The  location  x=  a^  of  these  perturbation  forces  is 

^ _ g-b  cos  20c 
~ (AL  ~ 2 cos  (z«<  + p) 


For  initial  angle  of  attack  equal  to  zero,  thla  becomes 


g-b 
2 cos 


It  should  be  noted  that  the  perturbation  may  be  caused  by  a change  of 
angle  of  attack  o(  , but  not  by  a change  of  flap  angle  s > since 
Is  zero.  ' 


dM/dp 


The  perturbation  center  of  pressure  a ^ Is  of  Importance  In 
unsteady  aerodynamic  problems  and  will  be  further  considered  in  the 
discussion  of  flutter. 
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Con5>aring  steady  and  perturbation  centers  of  pressure  aj^  and.  , 
it  is  seen  that  the  letter  is  in  general  well  upstream  of  the  former, 
a^  lying  in  the  vicinity  of  the  quarter  chord,  in  the  case  of  an 
elliptical  profile. 

Since  control  reversal  involves  the  change  in  pitching  moment 
about  the  elastic  axis  resulting  from  change  in  control  position,  it  is 
obvious  from  Equation  (8.1)  that  althovigh  a chaiage  in  ^ will  change 
the  lift,  it  will  not  effect  the  pitching  moment  about  the  midchord. 

For  elastic  axes  slightly  fonard  of  the  midchord  and  with  sufficient 
torsional  rigidity,  it  appears  that  the  problem  of  aileron  reversal 
should  not  cause  any  difficiilty. 


8.2  AERODYHAMIC  DERIVATIVES  FOR  FLUTTER  OF  AN  ELLIPTICAL  AIRFOIL 


Using  the  procedures  and  notations  of  Reference  (9),  the  complex- 
potential  associated  with  the  translation  and  rotation  without  circula- 
tion of  an  elliptic  cylinder  with  semi -axes  a and  b,  which  is  instan- 
taneously fixed  at  the  origin  of,  and  oriented  parallel  to  a coordinate 
system  x,y,  is 


JUT  = 

where 

a+  y 

A - 

t.io 

Co  = 

UO 
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Figure  8.2 


U, and  V,  are  the  velocity  components  in  the  positive  x and  y directions. 
We  define  the  complex  variables 


i ' 


$.11 


where 

Equation  (8.12)  is  the  transformation  equation  relating  the  elliptic 
contour  in  the  z plane  with  the  circular  contoiir  of  radius  A in  the  ^ 
plane,  ^ * A e^®. 

Since  this  is  a linear  theory,  we  can  separate  the  forces  and 
moments  on  the  elliptic  airfoil  into  non-circvilatory  and  circulatory 
parts  and  then  superpose  them.  The  former  is  obtained  exactly  by  use  of 
Lagrange's  Equation  for  moving  coordinates  (Reference  (9))  whereas  the 
latter  is  obtained  approximately  by  assming  proportionality  with 
certain  terms  of  the  Theodorsen  theory. 
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To  obtain  the  kinetic  energy  without  circulation  for  use  in 
Lagrange's  Equation  we  obtain  the  potential  function  and  the  stream 
function  on  the  contour  of  the  cylinder  from  Equation  (8.9) 


^onc 


U,  b cx)5  © -f  \(a  sin  0 


% 

— r2.c.b 
4- 


Sin  ZO 


or 


% 


one 


^ -Q-qV  rrm 

za.  J ' a'- 


SJ3b 


and 


one 


- - L/,b sin© 


-h  Va.  CO6  e cos,  20 

' • 4- 


e.14- 


The  kinetic  energy  is 

T = - 


^-V,a  + 


sCc.t  1 

<?  J 


The  upward  force  F and  counterclockwise  moment  L exerted  by  the  liquid 
on  the  moving  cylinder  is  given  by 


<^21  ^ ^ 
~ d±  91*  ^ 

©il  '9SL  % 


L 


where 


f - Mi 

A 

* SL  k 


?./7 
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We  obtain 


F = 


-Tr|5  - a.\av^)  + ^ (a'^v;  + b Ji q) J 


8.i<? 


In  order  to  observe  the  usual  sign  conventions  for  a cylinder 
moving  in  the  negative  dir??ction  relative  to  the  fluid,  vertical  deflec- 
tion positive  downvard,  lift  positive  downward,  angular  rotation  and 
moment  positive  nose  upward,  we  make  the  following  change  of  coordinates 


u,  - - u, 

V,  - - (li  S.ZO 

L - ~ '^neo 

F M = - P 

' J 
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We  then  obtain 


P„„  = - TTf  [a* ( 4 -bX'i:] 

-'^f[  ^ 


The  non- circulatory  force  Pjj^,  for  rotation  of  the  airfoil  about  an 
elastic  axis  x = a is  obtained  by  replacing  h by  h - aoc  , so  that 


= - Tr|>  a'  [(d  - a ^ 


8.22 


Similarly,  the  non-circulatoi^  moment  about  the  elastic  axis  x = a is 


8.2^ 


The  above  expressions  are  In  agreement  with  TheodorBen's  non- circulatory 
forces  and  moments  for  the  case  b =■  0,  a flat  plate. 

For  the  circulatory  forces  and  moments  due  to  the  wake  we  assume 
that  the  wake  vorticlty  is  transported  by  the  stream  so  that  It  possesses 
a constant  velocity  Uq  relative  to  the  airfoil.  We  further  assume  that 
the  induced  flow  about  the  airfoil  is  the  same  as  that  about  a flat 
plate  of  chord  2a., 

I 

We  now  introduce  a Kutta  condition  for  the  purpose  of  establishing 
the  strength  of  the  total  wake  vorticlty.  We  see  that  the  unsteady 
component  of  velocity  at  any  point  of  the  ellipse  can  be  expressed  in  the 
form 
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8.24 


I 


I - 


The  velocity  at  a point  x of  the  ellipse  due  to  the  total  vorticity  is 
of  the  form 


/ 


f,(5x)  U(5)de 


•ST 


8.25 


where  s is  the  element  of  length  of  the  trailing  streamline  measured 
from  the  separation  point.  By  equating  expressions  (8.24)  and  (8.25) 
the  point  of  the  ellipse  x = Xq  = a cos  p we  obtain  the  relationship 


at 


Q = ( h + + 


ITT 


<0 


g*  C ~ 2.  ) ~ O- 


8-2(0 


The  moment  about  the  reference  point  a due  to  the  vorticity  in  the 
wake  is  expressible  in  the  form 


= ZTT 


r <^0  00 


Q 


^^5 
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Mp  = ZTIj^iy^c'  6 (k,p)  Q S.27 

where  B(k, A ) is  the  ratio  of  the  two  infinite  Integrals  as  written 
above . 

Examining  the  non-cii^ulatory  moment.  Equation  (8.23)  we  see  that  it 
can  be  written  in  the  form 

Mhc  " •irj>c5L;.(h +U.OC)  + Tpa»[5  h - + -|i)  ] 


¥j>c*Ui[Q-  (a 


*.2  8 


Adding  Equations  (8.27)  and  (8.28)  we  see  that  the  total  moment  is 
written 


[ (of  Ctos'p-  1)-  a)Uo‘< 


Co 


a*  COS 


f 


-zh]  +2Trj’U.c:  [i  C„(^,  I,}  Q] 


where  al  will  be  seen  to  be  the  location  of  the  center  of  pressure,  and 
(^(  ^ ,k)  is  a new  function  of  ^ and  k to  replace  B(^  ,k). 

The  most  general  possible  form  of  1,  consistent  with  linear  theory, 
can  be  written 
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is 


i-  ^ L * * t-i.« 

But  Tq  to  "be  linear  requires  that 

Iq  = i.Q 

If  we  go  to  the  limit  (*)  — > O , we  find  that  steady 

- 21T  UoCjp,0)]  cl  «,30 

From  Equation  (8.2)  of  the  previous  section,  we  find  that  the  "steady" 
moment  due  to  small  perturbations  of  the  angle  of  attack  from  an 
initially  zero  value  is 


h 


so  that 

Cjto  ^ CL  CO^f> 


i ^ 

^ \ 

2.  COS  (3/ 


Substituting  the  value  of  Jlo  obtained  from  Equation  (8. 31)  into  (8.29), 
we  obtain 


[,  ^4  . £a.(cos fl— t)-a 


+ 2Tr 


g,32 
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‘Ihe  circulatory  lift  is  obtained  in  a similar  manner.  Thus 


where 


Cp(-p,k)  = 


Then,  adding  Equations  (8.22)  and  (8.33) 

p=  [(A'+-ai  ) + 

- iTT^l^o  ^ ( jS,  k)  Q. 


»-34- 


When  g and  b are  zero  Equations  (8.32)  and  (8.35)  reduce  to  Theodor- 
sen’s  expressions  for  the  flat  plate. 

From  Equation  (8.32)  we  see  that  the  first  term  in  the  bracket  is 
conservative  and  does  not  contribute  to  energy  exchange  whereas  the 
second  and  third  terms  are  unconservative  and  contribute  to  the 
mechanism  of  flutter.  The  second  term  vanishes  when 


a - 


S.3C, 


and  the  third  term,  when 


a.  - 


I0  - a. 

2 C06  6 


8,3? 
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For  the  flat  plate,  a - 0 and  = a.  Equation  (8.36)  corresponds  to 
the  three-quarter  ch(&d  and  Equation  (8.37)  to  the  quarter  chord.  In 
the  case  of  a thin  ellipse,  and  ^ small,  there  should  not  he  much 
change  in  these  locations.  Since  the  criterion  for  precluding  flutter 
is  associated  with  the  qviarter  chord  location,  it  appears  that  the  same 
will  he  true  for  the  thin  ellipse. 

In  the  case  of  orthodox  airfoils  incorporating  control  svirfaces 
such  as  ailerons,  control  reversal  inlts  sloiplest  aspect  may  he 
thought  of  as  arising  hasically  from  the  added  load  due  to  control 
surface  movement  acting  in  the  vicinity  of  the  midchord,  while  added 
load  due  to  incidence  changes  acts  near  the  quarter  chord.  Because  of 
the  latter  centroid  of  loading  it  is  usual  practice  to  locate  the 
elastic  axis  of  an  orthodox  system  in  the  vicinity  of  the  quarter  chord. 
We  see  then  that  the  use  of  a flap  in  the  usual  airfoil  results  in  the 
loss  of  control  surface  effectiveness.  Any  helicopter  hlade,  incorpora- 
ting the  usual  design  would  experience  hoth  loss  of  control  effective- 
ness and  undesirable  vihration  in  forward  flight.  Since  the  Thwaites 
airfoil  incorporates  a small  chord  flap  located  in  a region  of  conpara- 
tively  small  loadings  and  induces  lift  hy  change  in  angle  fS  so  that 
there  results  no  shift  of  center  of  pressure,  the  control  forces  and 
steady  forces  act  at  the  same  location,  thereby  allowing  for  increases 
in  lift  or  control  forces  without  loss  of  control  effectiveness.  To 
accomplish  this  it  is  necessary  only  to  locate  the  elastic  axis  in  the 
vicinity  of  the  midchord.  An  alternative  and  possibly  more  practical 
means  of  varying  lift  forces  on  the  Thwaites  airfoil  would  he  to  incor- 
porate a flap  which  is  hinged  at  a fixed  point  on  the  ellipse.  The 
cyclic  variation  of  force  could  he  obtained  hy  oscillating  the  flap 
about  its  root  hinge,  relative  to  a static  position.  Although  this 
would  require  a somewhat  larger  flap  chord,  it  is  believed  that  the  re- 
sulting load  would,  as  in  the  case  of  an  orthodox  flap,  still  act  in  the 
vicinity  of  the  midchord. 

Although,  according  to  previous  arguments  the  Thwaites  system 
appears  to  have  good  characteristics  from  the  standpoint  of  reversal  and 
divergence,  its  characteristics  with  reference  to  flutter  are  not  as 
clear  cut.  Assuming  that  an  elliptical  airfoil  is  used,  it  appears  that 
since  the  static  center  of  pressure  is  at  the  midchord  and  the  perturba- 
tion center  of  pressure,  associated  with  flutter  is  in  the  vicinity  of 
the  quarter  chord,  that  an  obvious  difficulty  arises.  It  is  believed 
that  by  incorporating  large  torsional  rigidity  while  maintaining  low 
flexural  rigidity,  a sufficiently  high  flutter  speed  could  be  obtained. 
However,  it  should  be  noted  that  the  elliptical  profile  has  been  the 
basis  of  the  present  analysis  only  for  piirposes  of  convenience.  It  is 
believed  that  other  profiles  of  an  oval  shape,  while  maintaining -the 
features  described  in  this  report  would  allow  for  the  two  centers  of 
pressure  to  be  closer  together. 
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Figtire  2.4  Velocity  Distribution,  Cj^  = 2 
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Figure  2.6  Velocity  Gradient,  Cj^  =i  1 
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Figure  2.7  Velocity  Gradient,  Cj^-  2 
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Figure  2.8  Velocity  Gradient,  Cj^=  4 
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Figure  3*2  Momentian  Thickness 
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Figure  3*3  Momentum  Thickness 
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Figure  3*5  Suction  Velocity,  Cj^  = 2 
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Figure  3*6  Suction  Velocity,  Cj^  - 4 
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Figvire  5*1  SpanwlBe  Variation  of  Blade  Angle 
and  Inflow  Angle  for  Plaseckl  Rotor 
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